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Abstract 

Today graphene nanoribbons and other graphene-based nanostructures can be synthesized with 
atomic precision. But while investigations have concentrated on straight graphene ribbons of fixed 
crystal orientation, ribbons with intrinsic curvature have remained mainly unexplored. Here, we 
investigate electronic transport in intrinsically curved graphene nanoribbons coupled to straight 
leads, using two computational approaches. Stationary approach shows that while the straight 
leads govern the conductance gap, the presence of curvature blurs the gap and reduces on-off ratio. 
An advanced time-dependent approach shows that behind the fagade of calm stationary transport 
the currents run violently: curvature triggers temporally and spatially focused electric currents, 
to the extent that for short durations single carbon-carbon bonds carry currents far exceeding the 
steady-state currents in the entire ribbons. Recognizing this focusing is pivotal for a robust design 
of graphene sensors and circuitries. 
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I. INTRODUCTION 


Graphene is a nanomaterial that can be tailored extensively for device fabrication pur¬ 


poses. It has already been used in field-effect transistors [U], non-volatile memory elements 
logical switches^, among other electronic components^]. A central part in the mak¬ 
ing of graphene leads, components, and circuitries is the patterning. Today experiments 


s, 


can routinely fabricate graphene nanori 
curvatures—even with atomic precision. 


Dons (GNRs) in a variety of shapes, widths, and 
6] For this reason especially straight GNRs have 
been investigated much both theoretically and experimentally. 5|, 7-ll| However, in addition 
to the straight GNRs, patterned graphene nanoribbons can also contain intrinsic curvature. 
Curvature has been demonstrated as semicircular graphene ribbons grown on a tern- 


plated silicon carbide substrate[12], as kinked and folded graphene with well-defined kink 
angles ( 9 ], and as nanoribbons sculptured by meniscus-mask lithography with customized 
edge curvatures 13]. In addition, curvilinear graphene microcircuits or stamps have been im¬ 


printed directly on graphene oxide films |l4|. But even though experiments have frequently 
realized curved structures, their transport properties have been remained unexplored. The 
stationary conductance simulations by Wurm et al. Q, Yin et al. 16], and Qiu et al. HQ 
are among the few investigations on curvature effects. We think that the additional degrees 
of freedom created by curvature deserve far more attention than previously invested. Ex¬ 
tending the views beyond straight graphene ribbons would help to improve the design of 
devices and nanocircuits. 

In this work we investigate electronic transport through curved graphene nanoribbons 
(CGNRs) using two theoretical approaches. First, we model CGNRs using density-functional 
tight-binding and calculate their transport properties by the regular Landauer approach [Q. 
These calculations are performed for a systematic set of CGNRs with different ribbon widths, 
edge curvatures, and lead couplings. Second, with the same set of CGNRs, we simulate trans¬ 
port using an advanced dynamic approach 20|, [21], a recent generalization of the Landauer- 
Biittiker formula for time-dependent transport 22f]. This approach provides a transparent 
analysis of the spatial and temporal dependence of electric currents, including their tran¬ 
sient dynamics. The analysis reveals that the stationary approach overlooks the complex 
processes that occur during hundreds of femtoseconds after switching on the bias voltage: 
in the time-dependent picture the currents in CGNRs show strong focusing both spatially 
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and temporally. 


II. CURVED GRAPHENE NANORIBBONS IN TRANSPORT CALCULATIONS 


To create models for the CGNRs, graphene disks were first cut out of pristine graphene by 
setting the origin at a hollow site, removing atoms below a minimum radius Ri and above a 
maximum radius R 01 and then removing the resulting singly-coordinated edge atoms. These 
disks were then bisected to yield 60°, 90°, or 180° arcs, to constitute the central scattering 
part for the transport calculation. Finally, to create the lead electrodes, the arcs were 
coupled to straight, semi-infinite graphene nanoribbons. The leads were of armchair type 
(AGNR) or zigzag type (ZGNR), depending on the bisection angle 23]. The 60° and 180° 
arcs were coupled to two AGNRs and the 90° arcs to AGNR on one side and to ZGNR on 
the other side. 

The three-part systems were saturated by hydrogen to remove the dangling bonds. Note 
that the systems have no in-plane stresses. We follow the convention to call the three parts 
the source electrode (S), the central conducting device (C), and the drain electrode (D) 

(Fig. [[]). 

We focus our discussion on a representative set of four CGNR samples (Fig. [[]). The 
samples include two 60° arcs of different curvature (samples V a and Vb, Fig. [T^). They both 
have 8-AGNR legs but in V a the curved parts are shorter and only the outer edge contains 
zigzag section; in Vb both inner and outer edges contain zigzag sections. The sample U 
contains a 180° arc and the sample L a 90° arc (Figs. [Do and c). These four samples are 
representative to illustrate the effects in transport properties caused by curvature. However, 
while the discussion focuses on these four samples, the general results are supported by cal¬ 
culations and detailed characterizations of a systematic set of different CGNRs, as presented 
in Supplementary Information (SI). 24] 


III. COMPUTATIONAL METHODS FOR STATIONARY AND TIME- 
DEPENDENT TRANSPORT 

To model the stationary transport, we used density-functional tight-binding (DFTB) 
method( 25 !- DFTB has been shown to capture the essential features of graphene-based 
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FIG. 1: (Color online) Curved graphene nanoribbons in transport calculations. Left panels: 
Schematics of curved ribbons attached to straight ribbons of the same width, forming a typical 
two-terminal transport device: source (S) and drain (D) electrodes connected via a central con¬ 
ducting device (C). Right panels: Atomic structures of the samples V a , Vb, U, and L. Hydrogen 
atoms are omitted from the figure for clarity. 


materials well 26T31I]. To prepare the samples, CGNRs were coupled to semi-infinite leads 
(lengths 88 A were sufficient), all edges were hydrogen-passivated to avoid finite-size 
effects, and the systems were relaxed. Self-consistent DFTB calculations then gave the 
T-point Hamiltonians for the curved parts, two principal layer units of source and drain 
electrodes, and their coupling matrices. These Hamiltonians were then used to construct 
the Green’s functions and self-energies for the three parts (L, C, and R) in order to calculate 
the stationary conductance by the Landauer formula. 19] 


To model the time-dependent transport, we used the advanced approach of Refs. 


20 


2l|. In this approach Kadanoff-Baym equations 32] are solved analytically to yield a time- 


dependent generalization of the Landauer-Biittiker formula 22]. Just like the stationary 
approach, also this approach assumes noninteracting electrons in a system composed of a 
central scattering part coupled to leads that are described within wide-band limit. Initially 
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all the system parts were in thermal equilibrium at a certain chemical potential p. At a 
given time the bias voltage Vsd = Vs — Vn was switched on symmetrically around /i. After 
the switch-on the reduced one-particle density matrix p{t) was evaluated directly from the 


closed-form expression with a time parameter t. 


211 ] Remarkably, the solution for p(t) appears 


in closed form and requires no time propagation. 

This approach provides a powerful tool for a transparent analysis of electron dynamics. 
Expressing the density matrix in a localized basis (of the system’s Hamiltonian) the diagonal 
elements yield local charge densities in individual atoms, whereas the off-diagonal elements 
yield currents in individual bonds. The approach thus provides a direct access to spatial 
and temporal information of electron dynamics. 

While the stationary transport was calculated using the full multi-orbital DFTB Hamil¬ 
tonian, the time-dependent transport was calculated by a single 7r-orbital tight-binding 
approximation with the hopping parameter 7 CC = —2.7 eV. Originally we began by using 
the full DFTB Hamiltonian also with the time-dependent approach, but it only added com¬ 
plexity to the analysis by giving noisy and complex current profiles within the transient 
regime. The 7r-orbital approximation was chosen because it simplified the analysis but still 


captured the central features of the time-dependent dynamics, even quantitatively 
(See also Sec. IV Al l 


33|, 


34] 


IV. TRANSPORT IN STATIONARY PICTURE 


In Fig. [2] we show the conductance of the four CGNR samples as a function of energy. 
The transport character of CGNR is dominated by the semiconductor or metallic nature 
of the AGNR or ZGNR leads, only with a slightly modified transport gap. Still, there 
are significant differences in the stationary transport of curved and straight ribbons, when 
compared directly (SI). 24!] Differences result from a pronounced electronic scattering that 
takes place at the interfaces between armchair and zigzag sections. Q 

The transport gap increases when the angle between the leads is increased. The gap is 
smallest for samples V a and Vb because the leads are symmetric and the central part creates 
relatively little scattering. In samples L and U the current has to turn a larger angle, which 
results in more scattering in the curved part and a larger transport gap. 


In straight ribbons the conduction channels open and close abruptly 24(, but in CGNRs 
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FIG. 2: (Color online) Stationary conductance of the four CGNRs of Fig. [0 


V a V U 



-0.31 eV 


FIG. 3: (Color online) Electron wave functions for LUMO (lowest unoccupied molecular orbital) 
and HOMO (highest occupied molecular orbital) in the CGNRs of Fig. Q] 

the channel energies are broadened by the asymmetry between outer and inner edges. Thus, 
one effect of curvature is to blur the transport gap. The blurring can be understood in a 
simple picture: an incident electron wavefunction propagating in the lead at a certain energy 
deflects in response to the curvature and creates an energy broadening. This broadening 
depends on the atomic details of the curved part. 

The stationary conductance can be understood via electron wavefunction visualizations 
at the T point. The highest occupied (HOMO) and lowest unoccupied (LUMO) molecular 
orbitals reveal both localized and delocalized characteristics (Fig. [3]). States relevant for 
transport show wavefunctions delocalized across the entire system. However, sample Vb 
contains zigzag sections at both inner and outer edges and shows strong HOMO and LUMO 
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localizations at the inner edge. The electrons in short zigzag segments are trapped by the 
neighboring semiconducting AGNR leads, leading to donor-like states near conduction bands 
and acceptor-like states near valence bands. Moreover, larger angle between the leads creates 
stronger localization of the frontier orbitals in the curved section, which is part of the origin 
for the observed large transport gap in L and U samples. 


V. TRANSPORT IN TIME-DEPENDENT PICTURE 


Transient behavior of currents in curved ribbons 


The stationary results suggest that the currents may run less smoothly through curved 
ribbons than through straight ribbons. To get a more complete microscopic picture of 
the current dynamics, we used the advanced time-dependent Landauer-Biittiker formula, 
described in Section m This time-dependent picture does not reduce full transport charac¬ 
teristics to a mere number, the energy-dependent conductance, but gives a more transparent 
access to both temporal and spatial dependence of the current. 

In the time-dependent approach, the bias voltage Rsd was switched on at t — 0 and 
the transient dynamics was calculated until a steady state was reached. The bias voltages 
ranged from Vs D = 0.2 eV to 2.0 eV. The coupling strengths were chosen so that the 
dissipation rate to the semi-infinite leads was 7 = 0.1 eV, for both source and drain. For our 
samples these conditions required saturation times up to 1 ns (see Fig. [7] and Ref. Q). The 
transient dynamics were first analyzed as currents through two bridges, B1 and B2, which 
were calculated as the sum over individual bond-currents (Fig. [4]). 

From the moment the bias is switched on, the current at B1 grows rapidly and starts to 
oscillate over hundreds of femtoseconds until the oscillations damp towards a steady state. 
To reach B 2 the wavecrest needs to travel ~ 30 A, which is seen as a short delay in the 
bridge current. This delay implies current velocity of ~ 10 A/fs, which equals the Fermi- 


velocity in graphene, as expected. 


351] The transients at both bridges are characterized by 


slow oscillations superimposed by fast oscillations associated with multiple intra-ribbon and 
ribbon-lead state transitions. The slowest oscillations originate from lead-to-lead reflections 
created by the charge density wave after switching on the bias. These slow oscillations are 
visible in the long-time plots of bridge currents (Fig. [7|and SI 241]). High bias reveals another 
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FIG. 4: (color online) Time-dependent currents in CGNRs. (a) Top figure shows the atomic 
structure of sample V a with bridges B1 and B2. The four panels show time-dependent currents 
through B1 (left) and B2 (right) at Fsd = 1-0 eV (middle) and at Vsd = 2.0 eV (bottom), (b) 
Same as panels a for sample Vb- Dashed lines mark the steady-state currents. 


type of slow oscillation that originates from multiple scatterings of the charge density waves 
trapped in the curved part. 


In straight GNRs the currents are regular and flow in the same direction 21], but in 


CGNRs they are more complex. The curved parts act as strong scattering centers for 
the electronic current and cause direction reversal, especially at low bias. Such backward 
currents are evident when looking at animations of bond current (see SI 24]). The dynamics 
of backward currents in CGNRs are governed by two length scales: lead-to-lead and lead- 
to-curved part distances. These distances create rich interference patterns upon multiple 
scatterings at the interfaces. The interference results in dense frequency spectra of the 
Fourier-transformed (FT) currents (Fig. [5]). 

It is instructive to investigate some transitions more closely. For instance, the time- 
dependent currents of V b at Vsd = 1.0 eV exhibit well-defined oscillations at 30 fs period, 
and the FT spectra shows an intense peak at the corresponding energy a; ~ 0.15 eV. This 
transition can be identified by looking at the local density of states (LDOS) at energies 
suggested by the spectral function. The local density of states at energies 0 eV and —0.15 eV 
show elongated edge states at inner and outer edges of the curved part (Fig. [6]). Thus, 
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FIG. 5: (color online) Fourier transforms of time-dependent B2 currents of Fig. |4) Samples are V a 
(left panels) and Vb (right panels) at the shown biases. 

the observed electronic excitations at cu ~ 0.15 eV corresponds to edge-state transitions 
taking place inside the CGNR. This oscillation becomes masked if Vgo is increased because 
additional electronic transitions start taking place within the enlarged bias window. 

Here we note that the single 7r-orbital model compares well with the full DFTB model. 
Both models yield frontier orbitals of similar edge-localized nature and similar level structure 
(Figs. Oandini), even if quantitative differences exist. Our central results are not affected by 
the choice of the model. 

In Fig. [7] we compare the transient currents in samples V a , L and U. Currents are calcu¬ 
lated at VgD = 1-0 eV through bridges in the middle of the curved parts. As a result, smaller 
angle CGNRs yield larger oscillation amplitudes in the transient regime. The sample V a 
with the shortest arc reaches the steady-state current fast, within ~ 300 fs. This is because 
the curved part is relatively short, with a zigzag section only on the outer edge, which means 
a strong coupling between the lead electrodes and a short equilibration time. The sample U, 
on the contrary, has a curved part that is longer and more structured, which means weaker 
coupling between the lead electrodes: current scatters more, creating smaller current fluc¬ 
tuations due to interference between scattered currents, and ultimately a longer equilibrium 
time. Sample L, with the same reasoning, has an equilibration time intermediate between 
samples V a and U. 

The spectral analysis for sample U shows an intense transition at cu ~ 0.14 eV, followed 
by several weaker transitions at higher energies (Fig. [8]). Sample L shows richer Fourier 
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FIG. 6: (color online) Density of states (DOS) of sample Vb- Shaded region shows the 1.0 eV bias 
window. Right insets show the local density of states (colormaps) for the selected DOS peaks at 0 
eV and —0.15 eV. 

spectrum. Consider, in particular, the triplet of intense transitions at 0.10 eV, 0.16 eV, and 
0.18 eV. These transitions can be tracked down to the triplet-like edge states at 0 eV that 
originate from three zigzag-sections in the structure (see Fig. |T]and spectral function in SI 
2dj). Therefore, this triplet of excitations, as well as the other triplets at higher energies, 
correspond to transitions between the edge states at 0 eV and states at other energies. 

The above analysis shows how FT provides some insight to the transitions during current 
transients. However, it misses all the temporal information. FT is unable to provide either 
the times or the durations for the activities of certain transitions. FT peak intensities give 
sort of averaged-out information about how long durations certain frequencies are present; 
a more accurate temporal analysis however requires an alternative analysis tool. 


B. Temporal and spatial characterization: 
Wavelet analysis 


The alternative analysis tool that gives better insight into the complex back-and-forth 
currents in CGNR is Wavelet decomposition 


36 


. In this method the time-domain signal 


is convolved by an element of a set of basis functions called wavelets. For our current results 
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FIG. 7: (color online) Left: Short-time behavior of currents through samples V a , U, L at Vgn = 1.0 
eV. Right: The corresponding long-time behaviors, showing current saturation. Currents were 
calculated through a bridge in the middle of the curved section. 



FIG. 8: Fourier transforms of the time-dependent currents in Fig. [7] for (a) sample U and (b) 
sample L at Vsd = 1-0 eV. 


we tested three types of wavelets: Haar[39j, Morlet|40i. 41], and Ricker 42|. [43], the last one 


being the second derivative of a Gaussian function. We will show only results using Ricker 
wavelet, as it provided the cleanest analysis. 

The continuous wavelet transform (CWT) of the current I(t) convolved by a given wavelet 
ip is defined as 

c(s ’ T) = 7 = sLj it)r ( i ^ 1 ) dt ’ 

where ip* is complex conjugate of the wavelet shifted by r and scaled by the dimensionless 
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FIG. 9: Wavelet coefficients (colormap) as a function of scale (s) and time (r) for currents through 
bridges B2 in (a) sample V a and (b) in sample Vb at Vsd = 1.0 eV. 


parameter s. This transformation maps the data into s x r space, where r is related to 
time and s can be related to frequencies via f s = f c /(sA), where A is the sampling period 
and f c is the center frequency of the wavelet. This frequency is calculated by associating a 
periodic function of frequency f c with the respective wavelet in such way that the function 
approximately delineates the wavelet form. The frequency is related to real time directly as 
t s = 1/f s - With fs as our time unit, we used Ricker wavelets with a center frequency of 0.25 
fs -1 and sampling of A = 0.16 fs. 

The advantage of wavelet analysis is that it enables us to distinguish between rapid and 
slow current fluctuations by choosing the scale s. Currents can thus be analyzed simulta¬ 
neously at different frequencies and at different times. By searching and analyzing intense 
wavelet coefficients we can identify instants and frequencies for eminent and possibly in¬ 
teresting events. We calculated CWT on the currents by using the Wavelet Toolbox in 


Mat lab 


441 ] and scales large enough to scan all the frequencies relevant to our samples. 


C. Spatial character of time-dependent currents 

The wavelet analysis of bridge currents in samples V a and Vb reveals frequency-dependent 
periodic structures of intense patterns (Fig. [9]). For sample V a the most intense wavelet 
components around s = 13-15 correspond to periods of t s ~ 9.0 fs and energy of 0.46 eV. This 
energy agrees quantitatively with the fourier peak at uj ~ 0.45 eV (Fig. [5]). Wavelet analysis 
also gives the “life time” of this mode: intensity peaks repeat every At ~ 14 — 20 fs and they 
dominate over other modes until 300 fs after which intensity starts to fade away. The modes 
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FIG. 10: Snapshots of bond currents (black arrows) and charge variations (colormap; units of 
electron charge) for samples V a (a and b) and Vb (c and d) at Vsd = 1-0 eV and at the instants 
shown in the panels. 


with scales s > 35 are the low-frequency modes with u < 0.065 eV that correspond to the 
background oscillations seen in the long run simulations [24]. Similarly, for sample Vb the 
most intense wavelet componets around s = 41 correspond to t s k, 27 fs and energies ~ 0.15 
eV, in agreement with the spectral analysis (Fig. [6]). Intensity peaks repeat every At ~ 55- 
60 fs and remain pronounced up to ~ 650 fs. The general current features in Vb are simpler 
than in V a and most wavelet componets retain their relative intensities during the course 
of time. The differences in these features can again be attributed to the longer curved part 
in Vb that more effectively delays the transmission of the current. The observed timescales 
provide a valuable information to help optimizing the detection of such transient responses; 
one could use a probing wave with the same period At to discern the characteristic modes. 

As discussed above, the wavelet analysis gives the possibility to choose interesting instants 
to further investigate not only currents through chosen bridges but bond current patterns 
within the entire samples. Wavelet analysis also enables identifying current “fingerprints” 
for each sample. Let us investigate few current pattern snapshots for samples V a and Vb, 
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with instants chosen so that one instant involves an intense wavelet amplitude and the other 
instant does not (Fig. HOjl . In both samples currents follow straight paths along the armchair 
leads and go through the inner zigzag edge in the curved part. Atoms directly connected to 
source and drain show the largest positive and negative charge density variations, and the 
spreading of the charge depends on the applied bias, on the length of the curved part, and 
on the electronic character of the AGNR lead. In sample V a the charge profile has more 
internal contrast, with density variations occurring within the entire sample, but with an 
extremely sharp transition from positive to negative density variations. This is because the 
leads are metallic 8-AGNRs that enable quick spreading of the charge density waves from 
the electrodes all the way to the short curved part. In sample V b the leads are also metallic 
but the curved part is longer, creating a more gradual transition from positive to negative 
charge variations and a smoother localization of the current along the inner edge. 

However, the striking feature in the snapshots discussed above is the strong localization 
or focusing of the current. We analyze the current in the sample V a by choosing the instant 
t = 12.6 fs based on the occurrence of the first visible peak in the wavelet amplitude. By 
creating a simple list of bond currents at that instant, we observe an extreme spatial focusing 
of the current (Fig. fTIT) : as few as two bonds carry instantaneous currents that are almost 
250 % of the overall steady-state current through the sample. At this instant also many other 
bonds carry currents larger than the steady-state current. Another instant with smaller 
wavelet amplitudes shows much weaker and less inhomogeneous currents. Therefore, the 
curvature seems to generate electric currents that show high temporal and spatial focusing. 
This is our central result. Such focusing is absent in straight GNRs, which show temporally 
and spatially far more homogeneous currents. 2l|] 

The current focusing is a generic feature common to all curved samples. For example, 
corresponding to the currents in samples U and L (Fig. [7]), Fig.[12]shows the wavelet analysis 
and Fig. [13] selected snapshots of two instants. The wavelet amplitudes and snapshots show 
current focusing, which is particularly evident in animations (see SI |24|). Sample U shows 
strong currents along the armchair leads and along the inner edges of curved section. Sample 
L has one lead zigzag and the other lead armchair, and this asymmetry becomes clearly vi¬ 
sible. Charge density variations occur along the sample irregularly, which happens because 
the zigzag lead supports more current paths than the armchair lead, at least for this type 
of CGNR. The difference in the propagation speeds in the different leads quickly create a 
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FIG. 11: Bond currents in sample V a at Vsd = 1-0 eV as a function of bond index (indexing wrt. 
increasing current). Plotted are the instantaneous bond currents at t = 12.6 fs (red circles), the 
bond currents at the steady-state (blue diamonds), and the overall steady-state current through 
the sample (dashed line). Arrow points to two singular bond currents that instantaneously carry 
almost 250 % of the overall steady-state current through the sample. 
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FIG. 12: Wavelet amplitudes (colormaps) as a function of scale (s) and time (r) for currents 
through the curved sections in (a) sample U, (b) sample L, and (c) sample V a with an impurity 
site, all at Vsd = 1-0 eV. 

non-uniform and irregular charge variations distributions across the sample. 

The wavelet analysis can also provide information about impurities. We repeated the 
calculations with sample V a , but this time with a model for an impurity atom adsorbed on 
top of a carbon at the bridge B2. The impurity is modeled by an on-site energy of —0.2 eV 
and a modified nearest-neighbor hopping parameter y, = 0 . 87 cc = —2.16 eV. This is a crude 
model for an impurity that causes weaker bonding to its neighboring carbon atoms. As 
a result of the impurity, the features of wavelet amplitudes get notably perturbed (compare 
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FIG. 13: Snapshots of the bond currents (black arrows) and the charge variations (colormap; units 
of electron charge) for samples U (a and b) and samples L (c and d) at Vsd = 1.0 eV and at times 
shown in the panels. 

Fig. [9ti with Fig. U2b). The fairly regular wavelet patterns of the pristine sample V a get 
disrupted, as becomes evident already within the first 100 fs. For instance, the second large- 
intensity region in s — r -space shifts towards larger scales, indicating that this particular 
mode becomes red-shifted. This shifting was confirmed by FT spectrum, which showed a 
peak at 0.20 eV for the pristine sample and at 0.18 eV for the impurity sample. The shifting 
occurs with scales s ~ 21-39, corresponding to periods between 13.8 fs and 25.7 fs. However, 
the most important effects of the impurity are the changes in current paths. Animations show 


how the impurity affects the charge distributions within the entire samp 


the current paths from outer to inner edge of the curved section (see SI 24]). In CGNRs the 


e and even alters 


electric currents can be extremely sensitive to single impurities: if impurities get adsorbed 
to sites with focused currents, the entire current patterns may have to rearrange. This can 
drastically change sample conductance. 

VI. DISCUSSION AND CONCLUSIONS 


Today experiments, by using sophisticated microscopy techniques, are able to investigate 
the dynamics of nanoscale materials at nanosecond 46] and at picosecond[47] time scales, 
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even if not yet quite at femtosecond time scales. However, even if direct measurements at 
the femtosecond scale should still remain out of reach, the events and dynamics at that scale 
nevertheless reflect behavior also in longer time scales. Our theoretical approach provides 
a transparent and intuitive perspective into the femtosecond dynamics and it also scales 
efficiently with respect to system size; it can be used to study experimentally relevant systems 
for phenomena with scales ranging from femtosecond dynamics to steady-state properties. In 
the future we will extend the method for other elements and other types of time-dependent 
perturbations. 

The simulations using this method show that curvature in graphene nanoribbons causes 
significant scatterings via impurity-like states, degrading the conducting channels in the 
vicinity of the transport gap. These scatterings cause strong temporal and spatial focusing, 


as seen in the visualizations enabled by wavelet analysis. The temporal focusing re 
transient phenomena, also observed in other studies and even in straight GNRs[21 


ated to 


48|,|49|, 

highlights the necessity of time-dependent approach to quantum transport. The curvature¬ 
generated spatial focusing is of fundamental importance and could be exploited in various 
applications, such as in chemical sensing and in the design of graphene circuitries in general. 

Scanning transmission electron microscopy (STEM) images of graphene constrictions 
exposed to adatoms have revealed that their sensitivity, as characterized by susceptible 


edge-states, can be superior to ’bulky’ graphene. 


501 ] Edge-state detectors, such as solid- 


state nanopore sensors applied to ultrafast DNA decoding, are currently under intense 


research. 


51 


53| . Decoding with single base resolution is challenging [54], and one of the 


related critical challenges is a sufficient spatial and temporal resolution. In this regard 
curved graphene structures might be suitable systems: currents at their edges are intense 
and highly focused. 
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